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Abstract 

Eiectromagnetic form factors of the transition ir + j v i r t. — > A\ are 
calculated by QCD sum rules technique with the description of the 
pion in terms of the set of wave functions of increasing twist. Obtained 
results are compared with standard QCD sum rule calculations. 
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Recently it was suggested to study pion form factor at not very large momen- 
tum transfers by light-cone QCD sum rules [jl] which combine the description 
of pion in terms of the set of wave functions of increasing twist with the tech- 
nique of QCD sum rules 0. This approach gives a possibility to calculate 
a contribution of so-called Feynman mechanism to the pion form factor, in 
which large momentum transfer selects the configuration, in which one par- 
ton carries almost all the momentum of the hadron. In that paper Q it was 
shown that at least up to the momentum transfers of order 10 GeV 2 this 
mechanism remains important and it is possible to describe pion form factor 
in this region by Feynman mechanism only. The second mechanism of large 
momentum transfer is the hard rescattering mechanism, in which large mo- 
mentum transfer selects configurations with a small transfers size in which 
the momentum fraction carried by interacting quark remains an average one. 
The hard rescattering mechanism involves a hard gluon exchange, and can 
be written in the factorized form [[|, ||, ||. 

There exists a number of arguments in favour of that at large but finite 
momentum transfer (Q 2 ~ 1 — lOGeV 2 ) Feynman mechanism is dominant. 
Using QCD sum rule approach [§] it was found ||,f7fl that pion form factor 
at Q 2 ~ 1 — 3GeV 2 is saturated by Feynman type contribution. However, 
this method can not be used for higher momentum transfer due to increas- 
ing contribution of higher operator product expansion terms at large Q 2 . 
Nevertheless in [§,|| using a concept of nonlocal condensates it was ob- 
tained indications that Feynman type contribution is dominate at least up 
to Q 2 ~ lOGeV 2 . At the same time an attempt to describe the data of pion 
form factor at Q 2 > 3GeV 2 by the contribution of hard scattering only leads 
to conclusion that the low energy pion wave function is very far from its 
asymptotic form [JlTJ] and there was proposed a model for pion wave func- 
tion which has a peculiar "humped" profile which corresponds that the most 
probable pion configuration is the case when one of the quark carry almost 
all pion momentum. But in |TJ| was pointed that the wave function of such 
type corresponds to respectively soft gluon exchange even at Q 2 ~ lOGeV 2 . 

Here we discuss Feynman type mechanism for n — ► A\ electromagnetic 



pion form factors. The first calculation of this amplitude was made in jlO 



by using operator product expansion for three-points correlator in vacuum. 



2 



1 



Let us consider the correlator: 

T^(p,q) = J e^d'x < 0\T{jl(x),f u m iO)}\n(k) >, (1) 

where ft = d'y^u and H m ' = ^uj u u — -^dj u d is the electromagnetic current, 
k is momentum of pion. This correlator was used in |l] to study pion form 
factor. 

Leading twist operator gives the following contribution to T^ v : 

T^u{p, q) = uf du- (hp 2 - q 2 )g^ v - 2(1 - u)p^p v 

Jo {1 — u)p A + uq z \2 

+ (1 - 2u)(p lJL q u + q^p u ) + 2uq fl q u ) (2) 

where (p n {u) is twist-2 pion wave function. We use the following definition 
for twist— 2 and —4 two-particle wave functions of pion: 

<0\d(0) llxl5 u(x)\ir(k) >=if v k lt I 1 e- iukx {ip{u)+x 2 gi {u) + 0{x*))du 

Jo 

+fn { Xft ~ kx~ x2 ) Jo e ' iUkX 92(u)du + ... (3) 

Here g\ and g 2 are twist-4 pion wave functions. 

In this paper we use different models for pion wave function. The first 
on is asymptotical wave function. In || it was shown that at asymptotically 
large Q 2 the pion wave function of leading twist has the following form: 

^ s \u) = Qu{\ - u) (4) 

The attempt to describe pion form factor at Q 2 > 3GeV 2 by hard rescattering 
mechanism only leads to a conclusion that the form of the wave function 
is much different from asymptotical one and it was suggested to use the 
following model for pion twist-2 wave function (see fTT||): 

(p { f' z \u, n ~ 500MeV) = 30u(l - u)(2u - l) 2 (5) 

And the third model pion wave function is 

<fl BF) (u) = 6u(l -u)(l + A 2 ^[5(2u - l) 2 - 1] 

+A 4 — [21(2« - l) 4 - 14(2u - l) 2 + lA (6) 
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which was proposed by Braun and Filyanov in at low normalization point 
H ~ 0.5Gey, with coefficients A 2 = | and A4 = 0.43 which is in agreement 
with QCD sum rules for first moments of pion wave function and provides 
that at a middle point u = 0.5 

^(0.5) ~ 1.2. (7) 

This value is with experimental values of various hadronic coupling constants 
calculated be light-cone sum rules approach (see ji~3|). 
It is easy to check that (0) satisfies to Ward Identities: 

qvF^v = -U(P- ?) At- (8) 

Twist-4 quark wave functions g\(u) and (fetV) give the following contri- 
bution into the correlator (HI): 



U [ du 
Jo 



4<7i (u) (i(p 2 -q 2 )g^ 



((1 — u)p 2 + uq 2 ) 2 
-2(1 - u)ptfp v + (p^q v + g^)(l - 2u) + 2uq fl q l/ ) 
-4g 2 (u)((l - u) 2 p^p u + m(1 - u)(p^q u + q^) + u 2 q^q v ) 
-AQ 2 {u){2{\ - «)p M p v - (1 - 2u){p tl q v + g^) - 2itfyg„)] 

2^) \ 

where ^ 2 (w) = - Jq g 2 (v)dv. 

There are also quark-gluon twist-4 operators: 



< 0|rf(0)7 A1 7 5 g s G' Q/3 (ra)w(x)|7r(A;) > 

7T J 



+ [k,(g, a - =) - k a (g,f3- /./ e^<^>Da*L(«0. (10) 

< 0|d(0)7 At ^ s G' Q/ 3(t;a;)M(x)|7r(A;) > 
= k^xp - x a k p )-U I e-^+^Da^^ai) 

71 J 

+ (h(g^ - ^) - **W ~ k -^)\ Uj e- ik < a ^Da^ ± ( ai ), (11) 



where G a/ 3 = |£ a) g M j,Cr MI/ , Dai = da\da2da^8 (1 — «i — «2 — 0:3). 
The contribution of these operators (|T0| , [TT|) is: 



U i du 



((1 — u)p 2 + uq 2 ) 2 



[(2(1 - u)p^p v - (1 - 2u)(p fjL q u + q^pu) 



where 



-2itfyg v )A(ii) - -g^(p 2 - q 2 )B(u) + 2(p M g„ - q^p u )C(u) 



A(u) 



r 1 -^ da 3 
da\ I 

ii Ju-otx «3 

u — a\ 



*|](Q!i) +2*j.(Oi) 



+ 1- 2- 



«3 



($i.(a i ) + 2$ ± (a i )) 



;i2) 



(13) 



B(u) 



u Ju-ai Q3 



( li — «i 

a 3 



<f>\\(cti) 



(14) 



C7(«) 



r 1 -"! a?a 3 
aot\ I 

Ju—a\ OJ3 



$ !.(<*) + 1 



2 )*i in, 

«3 



(15) 



There are also 4-particle twist-4 operators which were considered in ref. [|1J . 
Contributions of such operators are not considered here. 

A systematic study of the higher twist wave functions ( excluding four- 
particle operators) has been done in the paper [13| . The set of wave functions 
suggested in the paper includes contributions of operators with lowest con- 
formal spin and also the corrections corresponding to the operators with 
next-to-leading conformal spin, which numerical values were estimated by 
the QCD sum rule method. This set is (hereafter u — 1 — u): 



$11 (a,-) 



$±(ai) 



120e5 2 («i — a2)ai«2«3, 
-1205 2 a 1 a 2 a 3 (- + e(l - 3a 3 )^ , 



30<T(ai - a 2 )«3 ( - + 2e(l - 2a 
30<5 2 (1 - a 3 )aj (- + 2e(l - 2a 3 ) ) , 



(16) 
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9i(u) 
where 



^-S 2 u 2 u 2 + -S 2 (uu{2 + 13mm) + 10m 3 (2 - 3m + \u 2 ) ln(u) 

+ 10m 3 (2-3m + -M 2 )ln(M) 
5 
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5 2 uu{u — u), 



(17) 



5 2 ~ 0.2GeU 2 , e~0.5. 



(18) 



The value of 5 was determined in [15 



Using this set of the wave function we obtain the following expressions 
for A(u), B(u) and C(u): 



A(u) 



B{u) 



C{u) 



<T (10m 2 (1 - u) 2 + e [-4m - 22m 2 + 52m 3 - 26m 4 - 4 ln(l - u 
+4u 3 (10 - 15m + 12m 2 ) In / ' 



u 



(19) 



5 2 (-10m 2 (1 - m) 2 + e [-4m - 22m 2 + 52m 3 - 26m 4 - 4 ln(l - u) 
+4M 3 (10-15M + 12M 2 )ln(^--^ ) , (20) 
205W(1-m) 2 (2m-1). (21) 



Let us consider hadron contribution into correlator T^ v . A\ -meson gives the 
following contribution into this correlator: 



m\ — p 2 



<0|jJ|i4i(p) >< A 1 (p)\jt m -\7r(k) > 



m\ — p 2 



g, v m 2 A G,{Q 2 )-V,Vu% (g 1 (Q 2 ) - (l - G 2 (Q 2 )) 

m A \ \ m A j j 



Ql 1 f G 1 (Q 2 ) G 2 {Q 2 ) 



m. 



m A/ 
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-2q„p v (g 1 (Q 2 ) + ^G 2 (Q 2 ) 

\ m A 

+q»q» (g,{Q 2 ) - (l - ?Q G 2 (Q 2 )) 



where 



< 0\jl\Ax(p) >= ie^m A f A , 



(22) 



(23) 



< A^Hk) >= \(2p-q,q)g uX 



(2 j9 -g) i ,g A )G 1 (Q 2 ) 



■((2p-?,g)^-g 2 (2p-g),)g A G 2 (Q 2 



(24) 



= 1.26GeV is Ax-meson mass, e\ is a polarization of the meson, f A = 0.2 

U T TY7V7Y 



which was determined in 16 from a fit to the absolute rate for r 



decay. Notice that in paper |L0 it was used another set for mass and A\- 

2 

coupling constant - ~ 6.0 and m A = I.VSGeV. This parameters were 
determined in Jl7|] by QCD sum rule, f A = a/2^ 1 . 

The heavier spin-1 states contributions have the same form as (]22]) . 

In the limit of massless quarks due to the conservation of axial current 
jp, the only massless state of spin-0 (pion) gives nonvanishing contribution 
into the correlator (]!]). This contribution is 



p- 



<0\fMp) ><n(p)\jt m -\7r{k) >-- 



%^(g 2 )(2 P - ? ) t 



(25) 



Pion form factor F n (Q^ was studied in 



To study G\ form factor we consider correlation function fi{p 2 ,q 2 )' 

T^u(p, q) = g^ufiip 2 , q 2 ) + ■■■(other tensors). (26) 

Then according to (|i~4"D and (p5|) it is possible to write the following disper- 
sional relation for the structure g^ v : 

f ATT&G^Q 2 ) + 1 /■« p(£, Q 2 )d.S 

m A — p 2 tc J s p 2 



/i(p2 g2) = + ir ^!M£ + (suMon t erms) , (27) 

mi — ?r 7r j so r> — s 
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where p(s,Q 2 ) is the spectral density of higher state (spin-1) contribution 
into the correlator (0), sq > m\ and pion contribution is absent. 

To suppress higher states contribution and to kill subtraction terms in 
we use so called Borel operator B which was suggested in fl2|: 



B P 2f(P 2 



(p2)(n+l) / 


d 


n\ [ 


dP 2 


p2 






p2 


T = n ' 





f(p 2 ) = f B (t), 



-p 



n — > oo. 



(28) 



Applying this operator to ( |2"7| ) we obtain 

1 roo 



f?(t, Q 2 ) = j A m\G x (Q 2 )e-^ + - / e~p{s, Q 2 )ds. 

71 Jso 



(29) 



Notice, that higher states contribution is suppressed by a factor e s ^ < 

g— so/* 

From other side, the leading twist-2 wave function gives the following 
contribution to f B : 

-l r>2, 



fB (— ) = k f 1 pMe-ygpfo 
2 Jo (1 — u) 2 

fn 



s + Q 2 



e ^s. 



(30) 



In (^) we represent the formula in the form of dispersional integral af- 
ter Borel transformation and ^f^fn ( s+ q2 j is equal to imaginary part of the 
correlator at p 2 = s. 

According to (|9|) and ([12]) twist-4 wave functions contributions are 



fi 



B{twist-A) _ fn f 1 ( Q 2 (4 : g 1 (u) - B(u)) 4pi(w) - B(u) 



2 Jo 



-B 



(— 

\s + Q 2 



_ uQ 

e (i-«)t du 



t{l-uf 

1 -k 

2 J« 

s + Q 2 _ J_ 
Q 2 t Q 2 



ds 



;i-u)» 



4 ^ (^f) 
s + Q 2 



s + Q 2 / 
e~~tds. 
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Due to asymptotic freedom in the limit s — > oo the imaginary part of the 
correlator tends to his perturbative value. So it is possible, like in usual sum 
rules, to estimate higher states contribution suggested that the imaginary 
part of the correlator ([I]) is equal to its calculable part (by operator product 
expansion) starting at s = Sq. Or by other words, here we use the standard 
concept of duality, which tell us that Ai— meson occupies the "region of 
duality" in the invariant mass s up to a certain threshold Sq ~ 3GeV 2 . 
Thus, to take into account higher states contribution we use the following 
limits in of integration (|30| , [3T1) : 



< u < — — -T7T or < s < s . (32) 
so + Q 2 

Using eqs. (|29| - p2D we obtain the following sum rule: 

rB(twist-2) AB(twist-A) 2 

G 1 (Q 2 ) = ± 7 ± 4 e^, (33) 

where we use a standard model for higher state contribution using the limits 
of integration fl32| ) in expressions for fi^ and f±^- 

Three models of the leading twist wave function (p n (u) were considered. 
In Fig.l it was shown Q 2 — dependence of G\ form factor which was calculated 
using eq.(|33D at t = 1.5GeV 2 and s = 3GeV 2 . 

Notice, that the dependence of Gi on the pion wave function is very weak 
at Q 2 ~ 2.5GeV 2 . The reason of this is that at Q 2 ~ s = 3GeV 2 we 
integrate a pion wave function in the region < u < 0.5. In the limit t — > oo 
this integral is equal to 0.5 and does not depend on the form of wave function. 
It is clear that using for respectively large t the dependence of sum rule on 
the pion wave function will be weak and at Q 2 ~ 2.5GeV 2 this dependence 
is compensated by a small changing of integration region over u. 

In this picture we show predictions which were obtained by Ioffe and 



Smilga in |T0[ from QCD sum rules for three-point correlator. Notice that 
there is a big disagreement between predictions of QCD sum rules for three- 
point correlator and the case of asymptotic wave function for 7r-meson. The 
best agreement with predictions of Ioffe and Smilga at Q 2 < 2GeV 2 we have 
in the case of BF-wave function. At Q 2 ~ 3GeV 2 the predictions of |T0 



two times smaller than the result obtained in this paper. Probably, this 
disagreement can be explained by large contribution of higher dimensions 
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operators which becomes large at this momentum transfer (see details in 
Hj). At Q 2 < 2GeV 2 disagreement between the case of BF-wave function 
and three-point sum rule is about 20% which is usual accuracy of QCD 
sum rules. In this picture we show the experimental value for Gi(0) which 



is determined from partial width of decay A\ — > rcy. According to |18 



T(A± — > 7T7) = 640 ± 2A6keV. Stability of the sum rule (|33|) is illustrated in 
Fig. 2. Notice that at Q 2 > 15GeV 2 continuum contribution becomes larger 
than 30%. Twist— 4 operator contribution is smaller than 10%. 



4 

To find G2 form factor, let us consider q^p v and g M g^ tensor structures of 
contribution of hadrons into the correlator (p. The only spin— 1 states give 
nonzero contribution into these structures. In the case q^Pu this contribution 
has the following form: 

Up\Q 2 ) = -Y.^-Jg?(Q 2 ) + %Gf{Q 2 )) , (34) 
~ mi — p \ mf J 

and in the case q^q u the contribution is: 

f,{v\ Q 2 ) = E {{ofm + ^\Q 2 )) - G«(Q 2 )) , (35) 

Y mi - p 2 \ \ mj J j 

where J2i is a sum over spin-1 resonances [A\ and higher states with J PC = 
l + ~, mi and fi are their masses and residues into axial current, and 
are their form factors. From other side, f a and fp were calculated in Section 
II. Thus using eqs. ( |34] , |35D we have the following sum rule for G<± form factor: 

- \Up 2 ,Q 2 ) - Up 2 ,Q 2 ) = E zJ^~i G f(Q 2 )- ( 36 ) 

The left side of this expression was found in Section II at p 2 ~ 1 — 3GeV 2 
and Q 2 > lGeV 2 Using Borel operator (see (p8[)) and using the model of 
continuum which was described in previous Section to take into account 
higher state contributions we obtain the following sum rule: 

fn f U0 - <2!« \ (\ + u)^(u) 

— e * / e 't 1 - 11 ) du < — 

J a Jo [ 1-u 
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(1 -u)H 



{{l + 2u){-2g x {u)-2Q 2 {u) 



+-A(u)) + 2u(l + u)g 2 (u) + C(u) 



(37) 



fn !4 

J A -/n 



so 



(3 S + g 2 )^ 



2(s + g 



2^2 



3g + Q 2 

qh [ s + g 2 

+ g 2 ) 




Q 2 dependence of G2 form factor for different wave functions are shown in 
Fig. 3. These results were obtained from sum rule (|38|) at t — 1.5GeV 2 and 
So = 3GeV 2 . It is interesting compare the results obtained with the first cal- 
culation of this formfactor which was made by Ioffe and Smilga They 
used QCD sum rules for three point correlator and their results are depicted 
in Fig.3 too. In the case of this form factor there is a good agreement between 
the case of asymptotic wave function and predictions which was obtained in 
[|Hj. We have a big disagreement between our results for CZ and BF wave 
functions at Q 2 ~ lGeV 2 . This disagreement can be explained by a large 
contribution of higher twist operators. At Q 2 > 2GeV 2 the contribution of 
twist-4 operators becomes smaller 20%. At lower Q 2 ~ lGeV 2 this contri- 
bution is more than 30%. Stablity of the sum rules is illustrated in Fig. 4. 
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In this paper we have used light-cone sum rules to calculate electromag- 
netic 7f — A\ form factors in the region of intermediate momentum transfers: 
lGeV 2 < Q 2 < 15GeV 2 . It was studied the dependence of the results ob- 
tained on the form of the pion wave function of leading twist. It was found 
that at Q 2 > 3GeV 2 the behaviour of the form factors in the case of asymp- 
totical wave function are very different from the cases of CZ and BF wave 
functions. It was shown that the results obtained in the cases of CZ and BF 
wave functions are not very different from each other at Q 2 ~ IQGeV 2 . 



11 



We compare the results obtained in this paper with the first calculations 
of these form factors which was made in framework of QCD sum rules for 
three point correlator by Ioffe and Smilga ]ID |. These two predictions are 



in agreement within 20-30% accuracy at Q 2 ~ 1 — 2GeV 2 . This accuracy is 
usual accuracy of QCD sum rules. It is shown that there are disagreements 
between these two approaches at Q 2 > 2GeV 2 and Q 2 < lGeV 2 which can 
be explained by large contribution of power corrections at high Q 2 in the 
case of three-points sum rule and by contribution of higher twist operators 
at respectively small Q 2 in the case of light-cone sum rules. 

We obtain that vector dominance does not work in the process. If we try 
to fit the form factors by p and p'-mesons then we find that //-meson contri- 
bution is dominant at Q 2 > lGeV 2 , and we can not use vector dominance 
to estimate p — A\ — 7r-coupling constant. 

It is shown that at Q 2 > 1 — 2GeV 2 form factor G\ increases with growth 
of Q 2 . This behaviour can be explained in framework of naive quark model, 
where Ai-meson is consist of two quarks in P-wave. 

Notice that it will be very interesting to measure G\ form factor at 
Q 2 ~ 2.bGeV 2 because here we have a very weak dependence on the form of 
pion wave function. The slope of the form factor gives additional information 
on the function. 

Acknowledgements. The author is grateful to V.M. Braun, N. Isgur, and 
A.V. Radyushkin for useful discussions. 



Appendix A 

In [[Uj it was shown that there are relations between wave functions of the 
same twist. Here one more relation for twist-4 wave function will be obtain. 
Notice, that a set of wave functions suggested in [|T^] is in agreement with 
this relation. It is a generalization of one of the relations obtained in |T~ 
Let us consider the following matrix element: 

— < 0| J(0)7 M 7 p 7 1 ,75M(x)|7r >= {g w q v + g vp q ii - g pu q p ) 

x jf due' iuqx [upM + ix 2 g 1 (u) - j^92(u)j 
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+ (gnpXu + gupXn - g^uXp) J q due luqx g 2 (u), 
and find first derivative of (|A1|) over x v 



WwM C due- iu ^{2 9l {u) + 2G 2 (u) - ug 2 (u)) 

J 

1 



) / due- mqx (2 gi (u) + 2Q 2 (u) + ug 2 (u)). 



o 



From other side, using equation of motion for massless quarks we have: 

Id 

< 0|(i(0)7 A1 7p7 ; ,75M(x)|7r > 



fndx 

l r 1 



1 f L 

— / vdv < 0|d(0)7^7 P 7 y 7 5 £ a #G^(t>a;)w(:r)|7r > 



1 rl 

= i— / vrfwa; a {(g w g v \ + gvpg^x - g^g P x) 
< 0|d(0)#G ai ,(ux)7A75w(x)|7r > +e^ pA < 0|d(0)^G Q;I ,7AM(x)|7r > 



= ^mp(^) f 1 vdvV ai e-^ +va ^ ($1,(0,) - 2$ ± («* 

« 

- /' vdvVate-^^ ($,|(oO + 2^(«0) • 

Using new variable u: 

u — ol\ 

Oi + WO3 = It 17 = 

eq. (|A5|) becomes 

^mp(^) / e~ iqxu duA{u) + z(^x p - x p g p ) f e~ i(!XU duB(r\), 

where 

^-"i do 



A(u)= dan — T (u-a 1 )($\\(a i )-2$ ± (a i )), 

JO Ju-at O3 v ' 

B{u)= da x \ —l^-ax) ($11 (a*) + 2# ± (c 

JO Ju-ai «q v 



'u—(Xl O3 
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Comparing eqs.( |A2| ) and ( |A7| ) we obtain: 
2g 1 (u) - 2Q 2 {u) - ug 2 {u) 

aa>\ I — y\ u ~ 

u Ju-ai a 3 



«1 



)(«||(Oi) 



2$ i fa,-: 



(A10) 



2g 1 (u) + 2Q 2 (u) + ug 2 (u) 

dan / —t(u- ai) ($|i(ai) +2* ± (a < )) . (All) 



The first relation ( |A10|) was obtained in ||13|| . The second equation (|A11| ) is 
a new one. Using ( |A10| ) and ( |A11| ) we obtain 



2&2 (u) + ug 2 {u) 



da 1 —f(u-a 1 )($ ± (a i ) + V ± (ai)). (A12) 

Ju-ax «3 



It is not difficult to check that the set of wave functions suggested in [13 



satisfies to this new relation (|A12 ). Notice that the set of wave functions 
was fixed by using two integral relations which are valid for any twist-4 
wave functions, and additional one which was obtained for asymptotic and 
preasymptotic twist-4 wave functions. Thus, this new integral relation is a 
generalization of the constraints obtained in [|13] for the case of asymptotic 
and preasymptotic wave functions. 

Let we write the last relation for twist-4 wave function which was obtained 
in O: 



92{U) 



da x ^ (2d> ± K) - d>| 

Ju—cei Ct3 v 



a, 



(A13) 



Equation ( |A12| ) can be useful for construction of a set of wave functions 
which are far from asymptotic ones. 
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